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1.0  INTRODUCTION 


This  study  is  to  determine  the  dependence  of  the  auroral  DMSP  electron  precipitation  events  upon  the 
geomagnetic  Kp  index.  Using  known  values  of  Kp,  this  study  determines  the  average  number  of 
occurrences  of  auroral  events  of  varying  intensity.  These  auroral  events  are  defined  by  an  electron  energy 
flux  (JEtot)  that  exceeds  specified  threshold  values. 

To  study  the  relationship  between  the  auroral  DMSP  electron  precipitation  events  and  the  geomagnetic 
Kp  index,  the  superposed  epoch  analysis  technique  was  implemented.  This  approach  was  used  by  Nagai 
in  determining  the  relationship  between  several  geomagnetic  indices  [Nagai,  1988].  Epoch  zero  is  the 
time  at  which  Kp  values  persisted  above  a  specified  threshold  value  for  a  specified  duration.  The  zero 
epochs  are  determined  from  the  available  data,  and  the  Kp  series  are  superposed  over  a  time  window 
ranging  from  one  day  before  the  zero  epoch  to  three  days  after.  The  superposed  Kp  series  are  then 
averaged  over  all  epochs.  The  same  procedure  is  used  to  obtain  the  corresponding  time  series  of  electron 
auroral  events  for  comparison  with  the  Kp  index. 

It  is  important  to  note  that  the  results  from  this  study  cannot  be  used  to  locate  these  auroral  events 
spatially.  Estimating  the  location  of  these  events  was  done  in  a  separate  study  on  the  same  DMSP  data 
[fiounar  and  Hardy,  1992].  The  present  study  can  be  used  to  estimate  the  auroral  activity  based  on  the 
geomagnetic  activity. 

In  the  next  section,  we  describe  the  approach  taken  in  generating  the  data  base  over  each  hemispheric 
pass.  We  then  discuss  the  Kp  index  and  the  number  of  events  as  a  function  of  the  Kp  threshold  and  Kp 
persistence  values.  We  end  this  report  with  some  suggestions  for  further  studies  that  can  be  very 
beneficial  in  understanding  the  auroral  activity  temporal  variations. 


2.0  DATA  BASE  GENEKAITON 


The  data  base  was  generated  by  processing  the  DMSP/F9  electron  precipitation  for  the  year  1990.  The 
satellite  is  in  a  circular  polar  orbit,  in  the  10:30  to  22:30  meridian,  and  at  about  840  km  altitude.  The 
electron  measurements  are  binned  into  four  electron  flux  ranges  and  the  magnetic  activity  Kp  value 
stored  for  each  hemispheric  pass.  The  electron  events  considered  are  those  with  average  energy 
exceeding  1  KeV.  The  four  energy  flux  bins  correspond  roughly  to  brightness  levels  of  2kR,  6kR,  20kR, 
and  60kR.  Note  that  the  electron  fluxes  were  recorded  in  KeV/(cm^-sec-ster)  and  are  converted  to  kilo- 
Rayleighs  (kR)  using  a  3969  A  spectral  line  emission  [Figures  12-19,  Chapter  12,  Jursa,  1985]. 

In  Figure  1,  the  Kp  values  observed  in  the  first  30  days  in  1990  and  the  corresponding  electron  events 
are  shown  as  a  function  of  time  in  days.  A  quick  examination  of  these  curves  reveals  a  definite 
correlation  between  the  geomagnetic  index  Kp  and  the  number  of  electron  events  that  are  observed  over 
several  hemispheric  passes.  A  complete  set  of  figures  of  the  Kp  values  and  the  number  of  precipitating 
electron  events  in  1^  is  provided  in  Appendw  A. 
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Figure  1.  Kp  index 


3.0  APPROACH 


For  a  given  Kp  threshold  value  (KTH)  and  a  specified  persistence  time  over  this  threshold  (KPT),  the 
oata  base  is  scanned  and  the  times  at  which  these  Kp  conditions  are  satisfied  are  determined.  These 
times  are  taken  as  zero  epochs.  That  is,  the  time  at  which  Kp  was  above  the  threshold  value  KTH  for 
KPT  hours  is  the  zero  epoch.  Once  an  epoch  is  located,  the  scanning  of  the  Kp  data  is  resumed  after 
Kp  values  drop  below  the  threshold  value  K  IH.  In  this  way,  the  instances  when  Kp  persists  above  KTH 
for  over  KPT  hours  will  only  be  considered  once.  The  Kp  series  are  superposed  over  a  time  period 
ranging  from  one  day  before  the  zero  epoch  to  three  days  after.  The  superposed  Kp  series  are  then 
averaged  by  the  number  of  epochs.  The  same  procedure  is  used  to  obtain  the  corresponding  time  series 
of  electron  auroral  events.  Note  that  the  samples  from  the  data  base  are  not  available  at  uniform 
intervals.  Each  sample  represents  one  hemispheric  pass  in  time;  which  spans  about  51  minutes.  This 
means  that  there  is  on  average  one  hemispheric  pass  per  hour.  In  other  words,  there  is  one  data  sample 
per  hour.  Sometimes  data  gaps  occur  so  that  no  measurements  are  available  for  a  few  hours.  In  this 
analysis,  the  original  data  samples  per  hemispheric  pass  were  interpolated  (when  gaps  occurred)  or 
averaged  (when  more  than  one  sample  is  available)  as  necessary  to  obtain  samples  at  every  1.2  hours. 
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4.0  DISCUSSION  OF  THE  RESULTS 


Some  of  the  results  are  summarized  in  the  Table  1.  The  columns  1  through  6  give  Kp  parameters  , 
which  are  the  Kp  persistence  in  hours  (PKT  in  column  1)  above  a  specified  Kp  threshold  value  (KTK 
in  column  2),  the  Kp  peak  value  (KPP  in  column  3)  and  its  time  of  occurrence  (KTP  in  column  4),  the 
rise  time  of  Kp  from  the  threshold  value  to  the  peak  value  (ART  in  column  5),  and  the  fall  time  from 
the  peak  value  to  the  threshold  value  (AFT  in  column  6).  The  remaining  columns  give  information  on 
the  maxima  of  the  number  of  occurrences  of  electron  events  and  their  time  of  occurrence  for  a  JEtot 
threshold  of  2  kR  (#Ev  and  time  in  columns  7  and  8,  respectively),  for  a  JEtot  threshold  of  6  kR 
(columns  9  and  10,  respectively),  for  a  JEtot  threshold  of  20  kR  (columns  11  and  12,  respectively),  for 
a  JEtot  threshold  of  60  kR  (columns  13  and  14,  respectively). 


Whenever  a  rise  or  a  fall  time  is  larger  than  the  time  window  available,  it  is  denoted  by  asterisks. 
Whenever  a  time  of  occurrence  is  unreliable  because  of  poor  statistics,  a  question  mark  is  appended  to 
its  value.  Note  that  a  negative  time  for  the  occurrence  of  a  peak  means  that  it  occurred  before  epoch 
zero. 


Table  1.  Kp  and  Electron  Events  Characteristics 
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All  the  figures  of  Kp  and  the  number  of  events  are  shown  with  respect  to  time  away  from  epoch  zero. 
These  values  are  shown  uniformly  in  time  at  every  1.2  hours  even  though  the  processed  DMSP  data  is 
based  on  hemispheric  passes  that  do  not  occur  uniformally. 

The  results  reveal  certain  properties:  (1)  The  peak  values  in  Kp  and  the  number  of  events  increase  with 
increasing  Kp  threshold  values.  (2)  The  spread  or  the  width  of  the  Kp  values  above  the  Kp  threshold 
value  increases  with  increasing  Kp  threshold  value.  (3)  The  spread  or  the  dispersion  of  the  number  of 
events  above  the  mean  increases  with  increasing  Kp  thresholds,  but  stays  unchanged  with  increasing  Kp 
persistence  time.  (4)  The  peak  of  the  number  of  events  tends  to  occur  at  the  same  time  or  a  few  hours 
later  than  the  Kp  peak.  (5)  The  peak  of  the  number  of  events  increases  sharply  when  the  Kp  thresholds 
are  increased  beyond  the  mean  Kp  value,  and  then  levels  off  for  Kp  thresholds  greater  than  the  mean 
Kp  value  plus  1.0. 


4.1  Kp  CHARACTERISTICS 


We  first  discuss  certain  characteristics  of  the  geomagnetic  index  Kp  as  the  Kp  threshold  and  persistence 
time  is  varied.  These  conclusions  are  drawn  from  Figure  2,  and  are  summarized  in  Table  1.  For  a  Kp 
threshold  of  3.0,  and  a  persistence  time  of  at  least  3  hours  (Figure  2a),  a  negative  bias  in  the  Kp  peak 
of  about  1.2  hours  is  observed.  This  bias  increases  in  magnitude  as  the  persistence  time  is  increased. 
When  the  persistence  time  is  12  hours,  this  bias  is  as  much  as  6.0  hours  (Figure  2d). 

The  Kp  peak  values  inaease  with  inaeasing  Kp  thresholds.  For  a  Kp  threshold  of  4.0,  the  Kp  peak  is 
4.9  when  the  persistence  time  is  3.0,  and  is  5.4  when  the  persistence  time  is  increased  to  9  hours  (Figures 
2a  and  2c).  One  reason  for  this  increase  is  the  fact  that  the  more  disturbed  the  geomagnetic  activities, 
the  longer  they  persist. 

The  Kp  rise  time  from  the  Kp  threshold  to  the  Kp  peak  value  is  generally  smaller  for  higher  Kp 
thresholds.  For  a  persistence  time  of  3  hours  and  a  Kp  threshold  of  3.0,  the  rise  time  is  2.4  hours 
(Figure  2a).  When  the  threshold  is  increased  to  5.0,  the  rise  time  drops  to  1.2  hours.  Furthermore,  a 
shorter  Kp  rise  time  generally  leads  to  a  shorter  Kp  fall  time.  The  Kp  fall  time  is  the  time  it  takes  Kp 
to  drop  from  the  peak  value  to  the  threshold  value.  For  the  same  Kp  persistence  value,  the  fall  time  is 
16.8  hours  for  a  Kp  threshold  of  3.0,  and  is  8.4  hours  for  a  Kp  threshold  of  5.0.  Another  important 
feature  of  Kp  is  how  the  rise  time  and  fall  time  relate.  For  a  persistence  time  of  3  hours,  the  fall  time 
is  at  least  twice  as  long  as  the  rise  time.  This  ratio  drops  to  slightly  above  one  when  the  persistence  time 
is  increased. 


4.2  ELECTRON  EVENTS  CHARACTERISTICS 


The  characteristics  of  the  number  of  events  as  the  Kp  threshold  and  persistence  time  is  varied  are 
discussed  below.  These  conclusions  are  drawn  from  Table  1  and  Figures  3  through  7.  Consider  first 
events  with  JEtot  exceeding  2  kR.  The  peak  value  of  the  number  of  events  increases  as  the  Kp  threshold 
increasci  For  a  persistence  time  of  3  hours  (Figure  3a),  98  events  are  observed  for  a  Kp  threshold  of 
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Time  from  Epoch  Zero  (days) 

Figure  2(a).  Superposed  epochs  of  Kp  index  over  a  time  window  of  4  da)^  with  a  sampling  rate  of  1.2  hours, 
corresponding  to  a  different  Kp  threshold.  This  figure  is  for  a  Kp  persistence  time  of  3  hours. 
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persistence  time  of  6  hours. 
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Time  from  Epoch  Zero  (days) 

Figure  2(d).  Same  as  Figure  2(a),  but  for  a  Kp  persistence  time  of  12  hours. 
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Figure  3(a).  Superposed  epochs  of  the  number  of  events  per  hemispheric  pass  for  energy  fluxes  (JEtot)  greater  than  or  equal  to  2  kR  over  a 
time  window  of  4  days  with  a  sampling  rate  of  1.2  hours.  Seven  curves  are  plotted,  each  corresponding  to  a  different  Kp  threshold  (same  as 
Figure  2).  This  figure  is  for  a  Kp  persistence  time  of  3  hours. 
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Time  from  Epoch  Zero  (days) 

Figure  3(c).  Same  as  Figure  3(a),  but  for  a  Kp  persistence  time  of  9  hours. 


Persistence  of  Auroral  Activity 
. DMSP/F9  SateUite  - 


s'" 

Hit 

HA 

*-a 


m 

u 


V 


£  0 
0) 

gtSJ 
a  jS 

V  o 
O 

••H  D'CVI 

(5!^ 


u 

o 

«~l 

n 

V 

u 

00 

(w 

o 

ja 


to  r-  <H  <0 

00  j>  in  CM 


o  in  o  in  o  in  en 
CM  CM  CO  CO  in 


11  II  II  II  II  II  II 

A  A  A  A  A  A  A 


Pi  ft  0<  D<  Pi  P< 

H  ^  ^  M  ^  ^ 


CQ 

a 

o 

> 

H 

Id 

Stq 

flcvi 

dCO 

0 

S 


s:^uaAa 


13 


)  -0.5  0.0  0.5  1.0  1.5  2.0  2.5  3.0 

Time  from  Epoch  Zero  (days) 

Figure  3(d).  Same  as  Figure  3(a),  but  for  a  Kp  persistence  time  of  12  hours. 


3  as  opposed  to  102  events  for  a  threshold  of  S.  The  mean  value  of  the  number  of  events  observed 
over  all  hemispheric  passes  is  63  events.  The  peak  of  the  number  of  events  increases  sharply  when  the 
Kp  thresholds  are  increased  beyond  the  mean  Kp  value,  and  then  levels  off  for  Kp  thresholds  greater 
than  the  mean  Kp  value  plus  1.0.  For  a  persistence  time  of  3  hours  (Figure  3a)  and  for  Kp  threshold 
values  not  exceeding  Kp  mean  of  2.3,  the  peak  value  increased  by  about  11  events  each  time  the  Kp 
threshold  was  increased  by  ^  For  Kp  thresholds  greater  than  or  equal  to  3.0,  the  peak  variation 
increases  by  only  2  events  for  every  Kp  threshold  increase  of  ^ 

The  fall  time  for  the  number  of  events  is  at  least  twice  as  long  as  the  rise  time.  For  instance,  for  a 
persistence  value  of  3  hours  (Figure  3a),  and  a  Kp  threshold  of  5.0,  the  rise  time  is  13  hours  and  the  fall 
time  is  32  hours.  The  change  in  this  ratio  is  very  small  and  hence  negligible  as  the  persistence  time  is 
increased  (see  Figures  3b  through  3d). 

Figure  4  shows  results  of  events  with  JEtot  exceeding  6  IcR.  The  peak  value  of  the  number  of  events 
increases  as  the  Kp  threshold  increases.  For  a  persistence  of  3  hours  (Figure  4a),  38  events  are  observed 
for  a  Kp  threshold  of  3  as  opposed  to  47  events  for  a  Kp  threshold  of  5.  The  mean  value  of  the  number 
of  events  observed  over  all  hemispheric  passes  is  23  events. 

The  foil  time  for  the  number  of  events  is  at  least  twice  as  long  as  the  rise  time.  For  instance,  for  a 
persistence  value  of  3  hours  (Figure  4a),  and  a  Kp  threshold  of  5.0,  the  rise  time  is  12  hours  and  the  foil 
time  is  34  hours.  The  change  in  this  ratio  is  very  small  and  hence  negligible  as  the  persistence  time  is 
increased  (see  Figures  4b  through  4d).  The  rise  time  is  generally  about  ^  day  and  the  fall  time  varies 
between  1.0  and  1.5  days. 

The  results  of  events  with  JEtot  exceeding  20  kR  are  shown  in  Figure  5.  For  a  persistence  of  3  hours 
(Figure  5a),  10  events  are  observed  for  a  Kp  threshold  of  3  as  opposed  to  16  events  for  a  Kp  threshold 
of  5.  The  mean  value  of  the  number  of  events  observed  over  all  hemispheric  passes  is  about  7  events. 

Figure  6  shows  results  of  events  with  JEtot  exceeding  60  kR.  For  a  persistence  of  3  hours  (Figure  6a), 

4  events  are  observed  for  a  Kp  threshold  of  3  as  opposed  to  6  events  for  a  Kp  thres.'iold  of  5.  The  mean 
value  of  the  number  of  events  observed  over  all  hemispheric  passes  is  3  events. 


4.3  COMPARISON  OF  Kp  AND  ELECTRON  EVENTS 


Figure  7  shows  the  Kp  and  the  number  of  events  for  all  five  JEtot  thresholds.  The  Kp  threshold  is 
selected  at  4.0  and  Kp  values  stay  above  this  Kp  threshold  for  at  least  3  hours.  These  figures  show  how 
well  correlated  the  number  of  events  are  with  Kp.  At  about  %  day  before  epoch  zero,  Kp  rises  very 
rapidly  from  a  Kp  value  of  3.0  to  a  peak  value  near  epoch  zero,  then  the  Kp  values  drop  below  the  Kp 
value  of  3.0  in  about  1  day.  These  same  features  are  seen  in  the  number  of  events,  with  a  rapid  rise  then 
a  slow  foil  in  the  number  of  events.  The  rise  time  from  the  mean  value  to  the  peak  value  is  about  ^  day 
and  the  foil  time  is  about  1  day. 

To  measure  the  correlation  between  Kp  and  the  number  of  events,  the  normalized  correlation  coefficient 
is  computed  for  Kp  above  the  specified  threshold  value.  'Hie  correlation  coefficient  is  defined  as  the  ratio 
of  the  expected  value  of  the  product  of  Kp  and  the  number  of  events  to  the  product  of  their  standard 
deviations.  Assuming  stationarity  in  Kp  and  the  auroral  events,  the  time  average  is  used  instead  of 
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Time  from  Epoch  Zero  (days) 

Figure  4(a).  Same  as  Figure  3(a),  except  that  JEtot  ^  6  kR. 
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Figure  4(b).  Same  as  Figure  3(b),  except  that  JEtot  ^  6  kR. 
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Figure  4(c).  Same  as  Figure  3(c),  except  that  JEtot  ^  6  kR. 


0} 

hi 

V  & 

IT 

E-  2 
0 

S  N 


Sll 

W  A 

•-3 


« 


X 

P<  •*■> 


^  CO  »f3  C\}  ^ 

u 

u 

hi . 

r/Joi/50kftO»/50 

^cMMcocoT^-tenn 

^  II  II  II  II  II  II  II 
•  /\/\/\/\/\/\/\ 
dE  Ch  cii  ^  A  Ph 

p“< 


n 

<4J 

d 

o 

> 

» 

dj 

d 

o>co 

dw 

d(M 

d 


0 


s:^U0Aa 


18 


0.0  0.5  1.0  1.5  2.0  2.5 

from  Epoch  Zero  (days) 


Persistence  of  Auroral  Activity 
DMSP/F9  Satellite  —  Year  1990 


S'i^xiaAg 


19 


Figure  5(*).  Same  as  Figure  3(a),  except  that  JEtot  t  20  kR. 
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Time  from  Epoch  Zero  (days) 

Figure  5(b).  Same  as  Figure  3(b),  except  that  JEtot  >  20  kR. 
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Time  from  Epoch  Zero  (days) 

nwre  5(c).  Same  aa  Figure  3(c),  except  that  JEtot  i  20  kR. 
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s  Figure  3(d),  except  that  JEtot  ^  20  kR. 
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Figure  6(b).  Same  as  Figure  3(b),  except  that  JEtot  ^  60  kR. 
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Figure  6(d).  Same  as  Figure  3(d),  except  that  JEtot  ^  60  kR. 
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Figure  7(a).  Comparison  of  Kp  index  with  the  number  of  events  per  hemispheric  pass  for  energy  fluxes  (JEtot)  greater  than  or  equal  to  2  kR 
over  a  time  window  of  4  days  with  a  sampling  rate  of  1.2  hours.  The  solid  line  is  the  smoothed  number  of  events  which  are  the  3-point  averages 
of  the  samples  denoted  by  plus  signs  (+).  The  dashed  line  is  the  3-point  smoothed  Kp.  This  figure  is  for  a  Kp  threshold  of  4.0  and  a  Kp 
persistence  time  of  3  hours. 
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le  as  Figure  7(a),  but  for  an  energy  flux  threshold  of  6  kR. 
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Figure  7(d).  Same  as  Figure  7(a),  but  for  an  energy  flux  threshold  of  60  kR. 


ensemble  average  in  estimating  these  quantities. 


From  the  figures,  the  correlation  coefficient  is  about  0.98  for  events  with  JEtot  greater  than  or  equal  to 
2  kR  and  6  kR  (Figures  7a  and  7d),  and  drops  to  about  0.89  for  JEtot  >  20  kR  (Figure  7c).  The  results 
in  Figure  7d  are  not  reliable  because  of  insufficient  data  samples.  A  correlation  coefficient  of  0.7  to  1.0 
shows  that  the  two  time  series  are  highly  dependent.  Hence,  we  can  assert  with  certitude  that  the 
geomagnetic  activity  as  represented  by  the  Kp  index  and  the  auroral  activity  as  measured  by  the 
precipitating  electron  events  are  highly  correlated  indeed. 


5.0  CONCLUSION  AND  FUTURE  STUDIES 


The  observation  of  enhanced  Kp  is  generally  indicative  of  increased  number  of  bright  arcs,  llie  increase 
in  the  auroral  activity  of  precipitating  electrons  is  significant  for  only  a  limited  time  period.  During  this 
time  ,  a  high  correlation  between  enhanced  Kp  and  the  number  of  bright  arcs  can  be  achieved.  The 
waiting  time  to  increase  the  chance  of  observing  bright  auroral  arcs  after  enhanced  Kp  values  are 
observed  should  be  between  3  to  6  hours.  Further,  the  substantial  increase  in  the  number  of  observations 
of  intense  events  as  a  result  of  enhanced  Kp  lingers  a  maximum  of  36  to  48  hours.  After  this  time 
period,  this  number  drops  to  the  average  number  of  random  observations.  In  conclusion,  the  nearer  the 
observation  can  be  made  to  the  time  of  elevated  Kp,  the  better  the  chances  for  observing  arcs. 

This  study  also  shows  that  increasing  the  persistence  time  of  high  Kp  values  above  a  certain  threshold 
does  not  improve  the  likelihood  of  occurrence  of  intense  auroral  events.  That  is,  the  persistence  or 
lingering  of  the  enhanced  geomagnetic  activity  (as  seen  in  Kp  values)  does  not  mean  that  the  probability 
of  observing  bright  arcs  is  inaeased. 

Future  studies  that  can  provide  additional  understanding  of  the  auroral  activity  events  can  be  readily 
accomplished  using  the  existing  data  base.  Two  studies  can  be  implemented  using  the  present  data  base, 
which  would  certainly  reveal  important  features  in  the  electron  auroral  events  that  the  present  study  did 
not  cover. 

One  study  would  cany  out  a  superposed  epoch  analysis  on  the  number  of  electron  events  as  a  function 
of  specified  thresholds  of  the  number  of  events  and  how  long  they  persist  above  the  threshold  values. 
This  differs  from  the  present  study  in  which  the  epochs  are  obtained  based  on  the  Kp  threshold  values 
and  their  persistence  duration  above  these  thresholds. 

Another  study  would  cany  out  a  superposed  epoch  analysis  on  the  probability  of  occurrence  of  electron 
events  as  a  function  of  specified  thresholds  of  the  number  of  events  and  how  long  they  persist  above  the 
threshold  values.  For  each  hemispheric  pass,  the  occurrence  of  an  event  will  be  indicated  by  1  and  no 
occurrence  by  0.  This  study  would  then  yield  an  averaged  probability  distribution  of  occunence  of  events 
around  the  zero  epoch.  A  probability  distribution  of  the  auroral  electron  activity  as  a  function  of  time 
can  be  a  valuable  tool  in  estimating  the  persistence  of  auroral  events. 
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APPENDIX  A.  GEOM/VGNETIC  Kp  INDEX  AND  AURORAL  EVENTS 


The  geomagnetic  Kp  index  and  the  number  of  auroral  events  with  energy  fluxes  greater  or  equal  to  2kR 
observed  in  each  hemispheric  pass.  The  data  are  based  on  the  DMSP/F9  flux  measurements  in  1990. 
These  curves  are  shown  versus  day  number. 
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Figure  A-6.  Days  152  -  182. 
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Figure  A-11.  Days  302  -  332. 
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